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Neutron stars contain persistent, ordered magnetic fields that are the strongest known in the Universe. However, their 
magnetic fluxes are similar to those in magnetic A and B stars and white dwarfs, suggesting that flux conservation during 
gravitational collapse may play an important role in establishing the field, although it might also be modified substantially 
by early convection, differential rotation, and magnetic instabilities. The equilibrium field configuration, established within 
hours (at most) of the formation of the star, is likely to be roughly axisymmetric, involving both poloidal and toroidal 
components. The stable stratification of the neutron star matter (due to its radial composition gradient) probably plays a 
crucial role in holding this magnetic structure inside the star. The field can evolve on long time scales by processes that 
overcome the stable stratification, such as weak interactions changing the relative abundances and ambipolar diffusion 
of charged particles with respect to neutrons. These processes become more effective for stronger magnetic fields, thus 
naturally explaining the magnetic energy dissipation expected in magnetars, at the same time as the longer-lived, weaker 
fields in classical and millisecond pulsars. 
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1 Introduction 

The purpose of this talk is to present and discuss some of 
the physical processes that are likely to be relevant in de- 
termining the structure and evolution of magnetic fields in 
neutron stars, almost regardless of their (so far largely un- 
known) internal composition and state of matter. For this 
reason, I do not discuss the fascinating, exotic issues of 
quark matter, superfluidity, superconductivity, and the like, 
but emphasize the much more pedestrian concepts of stable 
stratification and non-ideal magnetohydrodynamics (MHD) 
processes such as ambipolar diffusion. A review with a very 
different focus has recently been given by Geppert (2006). 

2 Charged particles in neutron stars 

The name "neutron stars" incorrectly suggests stars com- 
posed exclusively of neutrons. However, additional parti- 
cles inside these stars play a crucial role. A neutron (n) in 
vacuum is known to decay by the weak interaction process 
n — > p + e + v e (beta decay) into a proton (p), an elec- 
tron (e), and an electron antineutrino (£> e ), with a half-life 
close to 15 minutes. This is impeded in very dense matter 
by the Pauli exclusion principle: If all the low-energy proton 
and electron states are already occupied, the only neutrons 
that can decay will be those with a sufficiently high energy. 
On the other hand, if there are many protons and electrons 



present, some of these will be energetic enough to combine 
into neutrons by inverse beta decay, p + e — » n + v e , where 
v e stands for an electron neutrino. In a neutron star, the neu- 
trons and protons will be confined by gravity, the electrons 
by the electrostatic potential of the protons (see, e. g., Rei- 
senegger et al. 2006), while neutrinos and antineutrinos are 
unbound and escape, contributing to the cooling of the star 
(e. g., Yakovlev et al. 2001). Direct and inverse beta decays 
will be in balance if the chemical potential^] of neutrons 
(Pn), protons (fi p ), and electrons (p e ) satisfy the relation 
M« — + Mei which forces the coexistence of a small frac- 
tion (few percent, but density-dependent) of charged parti- 
cles with a much larger number of neutrons (e. g., Shapiro 
& Teukolsky 1983). Additional particles (both charged and 
uncharged) can appear by other weak decay processes at 
densities higher than typical nuclear densities. 

In addition to stabilizing the neutrons, the charged parti- 
cles play two important roles regarding the magnetic fields 
and their evolution: 

- Being charged, these particles can generate electrical 
currents, which support potentially very strong magnetic 
fields. 

- Since the proton fraction Y = n p /n depends on den- 
sity (m stands for the number density of particle species 
i = n, p, e, and n = n n + n p is the total baryon density), 
neutron star matter is inhomogeneous, stabilizing it with 
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1 At zero temperature, these chemical potentials reduce to the respective 
Fermi energies. 
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respect to convective overturn (Pe thick 1992; Reisen- 
egger & Goldreich 1992; Reisenegger 2001a). As dis- 
cussed below, this is likely to have an important stabi- 
lizing effect on magnetic field configurations. 

3 Magnetic fields in neutron stars 

Many neutron stars are detected as pulsars, whose regu- 
lar pulsations in the radio, X-ray, and/or optical bands are 
produced by a strong magnetic field being turned around 
at the stellar rotation period P. These periods slowly in- 
crease in time, i. e., the neutron stars lose rotational energy, 
probably through magnetic coupling with their surround- 
ings. Modelling this coupling as electromagnetic radiation 
from a dipole rotating in vacuum, oriented orthogonally to 
the rotation axis, one can infer the surface magnetic field 
strength B oc \J PP, where P is the time-derivative of the 
rotation period (e. g., Shapiro & Teukolsky 1983). Inferred 
fields range from 10 8 G in millisecond pulsars up to 10 15 G 
in soft gamma-ray repeaters (SGRs; Duncan & Thompson 
1992; Kouveliotou et al. 1998; Woods et al. 1999), the latter 
being the strongest magnetic fields known in the Universe. 

In radio pulsars, the rotational energy loss can account 
for the whole observed energy output (relativistic particles 
and electromagnetic radiation) from these objects. For the 
strongly magnetized, but slowly rotating SGRs and anom- 
alous X-ray pulsars (AXPs), however, the observed X-ray 
luminosity is much larger than the rotational energy loss 
rate, so an additional source of energy is required, the most 
likely being the decay of their magnetic field (Thompson & 
Duncan 1996). This would make these objects be the only 
known magnetically powered stars, or "magnetars". An in- 
terest way of probing the strong magnetic fields inside these 
objects appear to be the quasi-periodic oscillations recently 
detected following two large flares of SGRs and interpreted 
as crustal shearing modes coupled to Alfven waves travel- 
ling through the stellar core (Levin 2007). 

In very old neutron stars, such as millisecond pulsars 
and low-mass X-ray binaries, the magnetic field is < 10 9 G, 
weaker than in young neutron stars, such as radio pulsars 
and high-mass X-ray binaries (~ 10 n ~ 14 G), suggesting 
that the magnetic field strength decays with time. However, 
this decay may be induced by accretion of matter from the 
binary companion (e. g., Payne & Melatos 2007 and ref- 
erences therein) rather than being a process intrinsic to the 
neutron star, so we do not consider it further. 

Magnetic field decay within the population of single ra- 
dio pulsars has also been suggested by some authors (Os- 
triker & Gunn 1969; Narayan & Ostriker 1990) but disputed 
by others (Bhattacharya et al. 1992; Faucher-Guiguere & 
Kaspi 2006), and does not seem to be well established. 

Aside from its intrinsic importance, the magnetic field 
may also play an indirect role in deforming the neutron star 
and in this way producing precession (Wasserman 2003), 
as appears to be observed in some pulsars, and gravitational 



waves, which might quickly reduce the rotation rate of new- 
born neutron stars (Cutler 2002). 

4 Origin of the magnetic field 

A natural hypothesis to explain the origin of the strong fields 
observed in neutron stars is the compression of the mag- 
netic flux already present in the progenitor stars. It led Wolt- 
jer (1964) to predict field strengths of 10 15 G, before any 
neutron stars had been identified observationally. Later, it 
was pointed out by many authors (e. g., Ruderman 1972; 
Reisenegger 2001b; Ferrario & Wickramasinghe 2005a,b, 
2006) that the distribution of magnetic fluxes is very sim- 
ilar in magnetic A and B stars, white dwarfs, and neutron 
stars, in this way providing support for the hypothesis of the 
magnetic fluxes being generated on or even before the main- 
sequence stage and then inherited by the compact remnants. 

On the other hand, Thompson & Duncan (1993) pointed 
out that newborn neutron stars are likely to combine vig- 
orous convection and differential rotation, making it likely 
that a dynamo process might operate in them. They pre- 
dicted fields up to 10 15 ~ 16 G in neutron stars with few- 
millisecond initial periods, and suggested that such fields 
could explain much of the phenomenology associated with 
SGRs and AXPs (Duncan & Thompson 1992; Thompson & 
Duncan 1995, 1996), some of which were later confirmed 
to spin down at a rate consistent with a strong dipole field 
(10 14 ~ 15 G; Kouveliotou et al. 1998; Woods et al. 1999). 

Of course, the two processes are not mutually exclusive. 
A strong field might be present in the collapsing star, but 
later be deformed and perhaps amplified by some combina- 
tion of convection, differential rotation, and magnetic insta- 
bilities (Tayler 1973; Spruit 2002). The relative importance 
of these ingredients depends on the initial field strength and 
rotation rate of the star. For both mechanisms, the field and 
its supporting currents are not likely to be confined to the 
solid crust of the star, but distributed in most of the stellar in- 
terior, which is mostly a fluid mixture of neutrons, protons, 
electrons, and other, more exotic particles. For this reason, 
the present discussion focuses on such an environment. 

5 Persistent, ordered field structures 

The magnetic fields of neutron stars, like those of upper 
main sequence stars and white dwarfs, appear to be ordered 
(with a roughly dipolar external configuration) and persis- 
tent (for much longer than a solar cycle, perhaps for the 
entire existence of these stars). As mentioned above, their 
magnetic flux distributions are similar, which also implies 
that their ratios of magnetic to gravitational energy (or mag- 
netic stress to fluid pressure) cover a similar range, up to 
~ 10~ 6 . Thus, even the strong fields of magnetars are fairly 
unimportant in terms of affecting the structure of the star. 
(The presence of superconducting protons might somewhat 
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increase this ratio, but only if their critical field is substan- 
tially larger than the actual field, which does not seem likely 
in the case of magnetars.) 

Another shared property of these stars is that a large part 
(if not all) of their interior is stably stratified; in the case of 
upper main-sequence envelopes and white dwarfs, because 
of a radially increasing entropy; in the case of neutron stars, 
because of a radial dependence in the fraction of protons, 
electrons, and possibly other particles. 

This makes it natural to search for magnetic field con- 
figurations that might be stable in stably stratified stars. An- 
alytic attempts in the past have failed, only yielding the gen- 
eral result that both purely toroidal fields and purely poloi- 
dal fields are unstable (Tayler 1973; Flowers & Ruderman 
1977), and the speculation that linked toroidal and poloidal 
fields might stabilize each other, yielding a stable equilib- 
rium (Prendergast 1956; Wright 1973). 

Recently, Braithwaite & Spruit (2004, 2006; see also 
Braithwaite & Nordlund 2006) carried out MHD simula- 
tions of stars with a stabilizing entropy gradient, which were 
initially given complex, "random" magnetic fields. These 
were observed to evolve on an Alfven-like timescale into 
a linked poloidal-toroidal configuration that persisted for a 
dissipative timescale. Thus, this configuration might be a 
good approximation to the actual field structures in upper- 
main sequence, white dwarf, and neutron stars. These would 
be stable in ideal MHD, but evolve over long times due to 
various dissipative processes present in each kind of star. 

In what follows, we will argue that the stable stratifica- 
tion of the stellar matter is crucial in allowing stable mag- 
netic equilibria. In neutron stars, the compositional stratifi- 
cation can be overcome by processes that either move dif- 
ferent particle species with respect to each other (ambipolar 
diffusion) or turn them into each other (beta decays). Both 
processes are enhanced by strong magnetic fields, yielding 
a natural explanation for the field decay in magnetars and 
its absence in other isolated neutron stars. 

5.1 Single flux tube as a toy model 

In order to understand the importance of stable stratifica- 
tion for the equilibria and evolution of the magnetic field, 
consider a horizontal magnetic flux tube with field strength 
B, fluid pressure P m , and fluid density p m embedded in a 
plane-parallel, unmagnetized fluid with pressure P ut> and 
density p out , stratified by a gravitational field. In order to 
be in mechanical equilibrium with its surroundings, it must 
satisfy P; n + B 2 /8tt = P out and p- m = p out . In a barotropic 
fluid, with a unique pressure-density relation p(P), these 
two relations cannot be satisfied simultaneously. If the first 
is satisfied, then P; n < P out and therefore p ln < p out , so 
the flux tube will buoyantly rise (Parker 1974). However, in 
neutron star matter, since the composition (here represented 
by the proton fraction Y) is inhomogeneous, the equation of 
state has the form p(P, Y). Thus, the flux tube can move to a 
position where a composition difference SY = Y in — Y out <~ 



'iB 2 / (%imp e ) compensates for the pressure difference, al- 
lowing the densities to be equal. So, the flux tube remains 
in equilibrium over time scales on which the matter can be 
modelled as a single fluid. 

Over longer time scales, two effects can erode the equi- 
librium described above. On the one hand, the matter inside 
the flux tube will be out of beta equilibrium, and this equi- 
librium can be restored by weak interaction processes such 
as the neutron beta decay or its inverse. On the other hand, 
particles can move in and out of the flux tube, in this way 
also equalizing the internal and external composition. On 
the characteristic time scale for these processes (which is 
very long, due to the strong degeneracy and particle colli- 
sions), the flux tube will move a distance over which the 
external composition changes by SY as defined above. This 
implies that the time scale for the flux tube to move a scale 
height is ~ Y/SY w 8irn e p e / (3B 2 ) times longer than 
the shorter of the beta equilibration and ambipolar diffusion 
times in the flux tube. 

In a more realistic geometry, one might consider a to- 
roidal flux tube in the equatorial plane of a spherical star. In 
this case, in addition to buoyancy, which tends to make it 
move radially outwards, the tube will be subject to a tension 
force, which exerts the opposite effect, namely to contract 
inward, towards the axis. At large radii, the former will be 
stronger, while the latter dominates at small radii. As for 
the plane-parallel case, this motion will actually happen in 
a barotropic star, but it will soon stop in one with a strong, 
stable stratification. In the latter, it can move only very little 
in the radial direction, so further motion is constrained to 
happen essentially on spherical surfaces. Thus, the flux tube 
will be unstable to contract on a spherical surface. 

5.2 Implications for stable equilibria 

As illustrated by the example of the single flux tube, the 
magnetic flux in the interior of a realistic, stably stratified 
neutron star can only move radially on long time scales 
controlled by weak interactions and ambipolar diffusion (as 
well as the Hall drift). In a combined toroidal-poloidal field 
configuration such as that found by Braithwaite & Spruit 
(2004, 2006; also Braithwaite & Nordlund 2006), the pres- 
ence of the poloidal field in the region around the magnetic 
axis impedes the toroidal field from contracting onto the lat- 
ter. At the same time, the poloidal field cannot deform sig- 
nificantly, because it is "tied together" by the toroidal field. 
Thus, it appears natural that such equilibria do exist in sta- 
bly stratified stars, but do not exist in barotropic ones. 

Additional insight can be reached through another, more 
formal argument: The combined pressure and gravitational 
force on a fluid element can be written as 

/ = -VSP + Spg, (1) 

where g is the gravitational acceleration, and SP and bp 
are small, Eulerian perturbations of pressure and density. 
If the whole configuration is axisymmetric, there will be 
no fluid force components in the azimuthal (0) direction, 
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thus the magnetic force in this direction must also vanish, 
j X B ■ <j> = 0. The "allowed" magnetic field configurations 
will therefore produce two independent force density com- 
ponents, in the poloidal plane, which must be cancelled by 
two components of the fluid forces of eq. (fl}. This is gen- 
erally possible in a stably stratified fluid, where SP and dp 
are independent functions of position. In a barotropic fluid, 
however, these two functions are proportional to each other, 
and the allowed equilibria will be much more restricted, if 
they exist at all. Thus, in a stably stratified fluid, there will 
be many equilibria, some of them likely stable, whereas in 
a barotropic fluid there will be few, if any. 



5.3 Implications for field evolution 

On the same time scales on which a single flux tube can 
move radially (controlled by weak interactions and ambipo- 
lar diffusion), the belt of twisted toroidal field lines in the 
Braithwaite-Spruit configurations can also move outwards, 
decreasing the magnetic energy stored in the star. However, 
it was shown above that the time scale for this to happen is 
very long. 

In fact, for the case in which weak interactions provide 
the dominant dissipation, the time scale to reach beta equi- 
librium is at best a few times shorter than the cooling time 
of the stat0 (Reisenegger 1995), unless the chemical im- 
balance r\ = fj, n — /j,p — (j, e ps B 2 /(8irn e ) is very large, 
1 77 1 S> kT. Since a substantial motion of the magnetic struc- 
ture takes - Y/5Y « 8im e n e /(3B 2 ) - (4 x 10 17 G/5) 2 
times the equilibration time, the field would never decay, 
were it not because the field decay itself could keep the star 
hot for a longer time (Thompson & Duncan 1996; Reiseneg- 
ger et al. 2005). A strong field will remain essentially con- 
stant until the temperature has decreased to kT w 1 77 1 / 5 - 5, 
at which point the weak interaction processes start having a 
net heating effect (Fernandez & Reisenegger 2005), con- 
verting the energy stored in the chemical imbalance into 
neutrinos and heat. Thus, the stellar interior remains essen- 
tially at this level until the magnetic field has decayed sub- 
stantially (Thompson & Duncan 1996). 

The ambipolar-diffusion-dominated limit is relevant at 
lower field strengths (the transition occurs near 10 16 G), 
and here again the heat input from the magnetic field can 
eventually offset the cooling and keep the star warm until 
the field has decayed. 

The time scales for these processes are difficult to esti- 
mate precisely (and may depend strongly on Cooper pairing 
effects, neutrino emission mechanisms, and the like), but es- 
timates based on Goldreich & Reisenegger (1992) appear to 
be roughly consistent with the time scales ~ 10 4 ~ 5 yr re- 
quired to account for the magnetar phenomenon. 



2 Both beta equilibration and cooling depend on weak interactions and 
have the same temperature-dependence, until the star is cool enough for 
photon cooling to become dominant. 



6 Crossing the solid crust 

Of course, neutron stars have a solid crust that must be 
crossed by the magnetic flux on its way out of the star. 
Its precise properties are still not clear, as it is a very un- 
usual solid, subject to a strong pressure, strongly stratified, 
and with a relatively weak shear modulus. Thus, it is not 
clear whether strong magnetic stresses would crack the crust 
in discrete events, suddenly releasing a substantial amount 
of energy and possibly accounting for the soft gamma-ray 
bursts (Thompson & Duncan 1996), or rather flow plas- 
tically, with a continuous release of energy (Jones 2003). 
If the magnetic forces are not strong enough to break and 
move the crustal solid, the flux will have to be transported 
by a combination of Ohmic diffusion and Hall drift (e. g., 
Jones 1988; Goldreich & Reisenegger 1992; Pons & Gep- 
pert 2007; Reisenegger et al. 2007). 

7 Conclusions 

The extremely strong magnetic fields found in neutron stars 
are nevertheless weak enough to be balanced by small per- 
turbations in the density, pressure, and chemical composi- 
tion of the stably stratified, degenerate, multi-species fluid 
found in their interior. The field likely originates from the 
fairly strong magnetic flux inherited from the progenitor 
star, possibly modified by a combination of convection, dif- 
ferential rotation, and magnetic instabilities acting during 
the short, protoneutron star stage following collapse. The 
equilibrium field set up during this stage likely involves 
linked toroidal and poloidal field components. It can decay 
through dissipative processes such as weak interactions and 
ambipolar diffusion, which change the chemical composi- 
tion of the matter, allowing it to move. These processes be- 
come particularly effective at high field strengths, possibly 
accounting for the energy release inside magnetars, which 
also keeps the stellar interior hot enough for the magnetic 
field to rearrange substantially. This forces changes as well 
in the crust of the star, which can be broken by strong fields, 
whereas weaker ones can evolve by a combination of Hall 
drift and resistive dissipation. The internal dissipation pro- 
cesses are ineffective for pulsar-strength fields, which can 
easily survive for a pulsar lifetime. None of these processes 
depends essentially on the exotic properties of a neutron 
star, such as Cooper pairing or quark matter, which can how- 
ever modify the time scales for these processes to occur. 
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